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A new point of view on the dynamics of the circumstellar nebula Ml-67 around the run-away Wolf-Rayet (WR) star WR 124 is 
presented. We simulated the outbursts of nebulae with different morphologies, to compare the results to the observed dynamical 
spectra of Ml-67. We found that it has been interacting with the surrounding ISM and has formed a bow shock due to its high 
velocity of about 180 km s^' relative to the local ISM. The star is about 1.3 parsec away from the front of this bow shock. 
The outbursts that are responsible for the nebula are assumed to be discrete outbursts that occurred inside this bow shock. The 
ejecta collide with this bow shock shortly after the outburst. After the collision, they are dragged away by the pressure of the 
ISM, along the surface of the bow shock. The bow shock is oriented in such way that we are looking from the rear into this 
paraboloid, almost along the main axis. Evidence for this is given firstly by the fact that the far hemisphere is much brighter 
than the near hemisphere, secondly by the fact that there is hardly any emission found with radial velocities higher than the 
star’s radial velocity, thirdly by the fact that the star looks to be in the centre of the nebula, as seen from Earth, and finally by 
the asymmetric overall velocity distribution of the nebula, which indicates higher radial velocities in the centre of the nebula, 
and lower velocities near the edges. We find evidence for at least two discrete outbursts that occurred inside this bow shock. Eor 
these outbursts, we find expansion velocities of Vexp ~ 150 km s“' and dynamical timescales of about 0.8 and 2 x lO'^ yr, which 
are typical values for LBV outbursts. We therefore conclude that Ml-67 originates from several outbursts that occurred inside 
the bow shock around WR 124, during an LBV phase that preceded the current WR phase of the star. 

Key words. Stars: circumstellar matter - Stars: individual: WR 124 - Stars: mass-loss - Stars: Wolf-Rayet - ISM: individual 
objects: Ml-67 - ISM: jets and outflows 


1. Introduction 


In this article, we describe our research on the dynamics of 
the Wolf-Rayet ring nebula Ml-67. Ml-67 is a bright nebula 
around the Wolf-Rayet (WR) star WR 124. The star has a high 
heliocentric velocity of almost 200 km s * and is also known 
as Merrill’s star ( |Vlerrill 1938) and 209 BAC. The star is clas¬ 


sified as a population I WN8 star (Bertola 1964) and is lo¬ 


cated in the constellation Sagittarius. Distance estimates vary 


from about 4.5 kpc ( 

Pismis & Recillas-Cruz 

1979 

) to 6.5 kpc 

(Nugis & Lamers 

2000|). The star has a terminal wind velocity 


is estimated to be about 20 Mq and its luminosity 6x10^ Lq 


(Nugis & Lamers 2000). 


The nebula Ml-67 around WR 124 shows a clumpy struc- 
ture, and most of the gas is concentrated in knots and filaments 


(Sirianni et al. 1998). An HST image of the nebula is displayed 


Send offprint requests to’. H.J.G.L.M. Lamers, e-mail: 
lainers@astro .uu.nl 


in Fig. ||. The nebula was first classified as an H ii region. After 
the discovery that the nebula has about the same radial velocity 
as WR 124, it was suggested that the nebula might be a plan¬ 
etary nebula (Minkowski 1946| ). However, the presence of a 
WR star and the N-enhancement and O-deficiency of the neb¬ 
ula suggest a WR ring nebula (Sirianni et al. 1998). The dis¬ 
tance estimates also point in the direction of an ejected nebula, 
so that Ml-67 is now generally accepted as a Wolf-Rayet ring 
nebula. 


Though WR ring nebulae are not necessarily ring-shaped, 
they often exhibit a structure of arcs or rings. This suggests 
that the nebulae may be created by discrete outburst events. 
It is generally thought that WR ring nebulae originate from a 
Luminous Blue Variable (LBV) stage of the central star, which 
is supposed to precede the WR phase. 

the LB Vs have circumstellar nebulae 
). The different nebulae are very similar 
in terms of physical properties. The expansion velocities are in 
the order of 50 to 100 km s ', their sizes about 1 parsec, and the 


More than half of 


(Nota & Clampin 1997 
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dynamical ages are in the order of 10"^ yr. The densities of the 
nebulae are generally found to be low (500 to 1000 cm“^) and 
the temperatures are in the range of 5000 to 10000 K (Lamers 
et al. 


It is still a point of debate whether the LBV outbursts occur 
during a Red Supergiant (RSG) or a Blue Supergiant (BSG) 


phase. The RSG scenario is proposed by Stothers & Chin 


(1993,1996), who suggest that the ejection of mass occurs only 
as a single event during a brief RSG phase. They explain the 
enhanced abundances of heavy elements of the ejecta by con¬ 
vective mixing in the RSG envelope (Stothers & Chin 1993, 


1996). 


On the other hand, Langer et al. (1994) suggest that after 


vent it from becoming a RSG (Langer et al. 1994; Lamers et al. 

zomh. 


of impacts of outbursts on the bow shock surface. In Sect. ^ we 
summarise the results and present the conclusions of this study. 

2. Observations 

For this research, we used the following three sets of observa¬ 
tional data. 


2.1. Long-slit spectra 

The first dataset we used, is velocity information from long slit 


spectra we obtained from A. Nota, published in [Sirianni et al 


the star has left the main sequence, it moves red-ward in the 
Hertzsprung-Russel diagram (HR diagram) and the expanding 
envelope becomes unstable, so that the star starts to develop 
extreme mass loss. This mass loss may be as high as 5 x 10“^ 
Mq yr“^ and is observed as the LBV ejecta. They explain the 
chemically enriched ejecta by rotation induced mixing. These 
eruptions therefore take place when the star is a BSG and pre¬ 


The research of the chemical composition of LB V-ejecta by 
Lamers et al. (2001) also indicates that the LBV eruptions oc¬ 
cur during a BSG phase. They suggest that the LBV outbursts 
are induced by the rapid, near-critical rotation of the star. In 
their scenario, the star is also being prevented from becoming 
a RSG by the mass loss. However, they point out that if a mas¬ 
sive, optically thick shell is being expelled from the star, it will 
cool as it expands and the physical conditions will temporarily 
be similar to that in the outer layers of a RSG, so that the for¬ 
mation of dust can also happen in this case. This mechanism 
explains the observed Humphreys-Davidson limit, that depicts 
the lack of RSGs with luminosities higher than 6 x 10^ Lq. 

The goal of our research is to disentangle the geometry and 
dynamics of the nebula Ml-67. In order to do so, we create 
different numerical models, of which the output is compared 
to available observations. Firstly, we present the observational 
data we use for this study in Sect. We will then discuss our 
models for freely expanding outbursts in Sect. The reason 
why we let these models expand freely is that the O-star that 
precedes a WR star blows a bubble of typically 3 0 pc in the ISM 
during its lifetime (Lamers & Cassinelli 1999). It will take an 

I _ j ^ I ' 

outburst of 100 km s more than 10 yr to cross this distance, 
so that it can indeed be considered to expand without any dis¬ 
turbance. We will find that no satisfying fit can be made, so 
that the assumption of a freely expanding outburst is wrong. 
We show that the cause for this is that the star has a high veloc¬ 
ity relative to the ISM. This causes a paraboloid-like bow shock 
instead of a more or less spherical bubble. The star is about 1 pc 
away from the front of this bow shock, so that an outburst of 
100 km s * needs only 10“* yr to cross it. Once it has done so, 
part of the outburst will collide with the bow shock, and will 
possibly be dragged away along the bow shock surface. We 
discuss the bow shock models for the case of a constant wind 
velocity in Sect. In Sect. |] we compare the theoretical bow 
shock models for a constant stellar wind with the observations 
and find a remarkable resemblance. Sect. ^ discusses the results 


( |1998|) . These data consist of 13 long slit spectra, taken with the 
ESO Multi Mode Instrument (EMMI) at the 3.5m NTT in La 
Silla. Each slit is positioned over the nebula in the east-west di¬ 
rection, at constant declination. The declinations of the slits lie 
between -30.82"and -i-24.69"relative to the star’s declination. 
In total, 413 good data points (right ascension, declination and 
radial velocity) were derived from the spectra, which formed 
the input for our study. These data points have been plotted in 

Fig.0. 

2.2. Fabry-Perot images 

The second datas e t are E abry-Perot (EP) images obtained from 
[Grosdidier et al. ( 1999 ). The EP images are shown in Pig. 0 
These 30 images were made in August 1996, using CPHT-SIS, 
with the etalon of the Universite Laval, Quebec, Canada. Each 
image was taken in Her, at a slightly different wavelength, so 
that it displays the emission at a certain radial velocity. These 
observations give a much more detailed view of the velocity 
distribution than the long slit spectra. Each image consists of 
100 X 100 pixels. Prom the darkest points, we used the coordi¬ 
nates (x, y, v) in the same way as the data points from the long 
slit spectra. A few thousand points were used. Note that the 
southern part of the nebula, more than approximately 20" south 
of the star, is missing (compare to Pig. 0 and see Pig. 0(b)) 
due to deteriorating seeing during the observations (Grosdidier 
et al. 1999). 


2.3. HST image 

Por reference , we also used the Hub ble Space Telescope im¬ 
age taken by Trosdidier et al. ( 1998 ), for example to identify 
structures found in the Eabry-Perot images. 

The image is a composite image of four WPPC2 images 
with a total exposure time of 10000 seconds, taken in March 
1997, using the narrow band P656N Ha filter. In the same way, 
four images of the same field were taken through the broad¬ 
band P675W R filter and combined to obtain a ’continuum’ 
image close to Ha. The continuum image was flux-scaled and 
subtracted from the first composite image to obtain the deep 
continuum-subtracted Ha image with the field stars removed, 
of which a negative version is shown in Pig. 0 

In this image different distinct arcs are visible. However, 
although these arcs are clearly seen locally, it is very difficult 
to find a global system of rings. The fact that the arcs are so 
clear, indicates that they are probably formed during different. 
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discrete outbursts, whereas the deficit of a global pattern may 
tell us that the history of the nebula is less straightforward than 
we might think a first. 


3. Models of freely expanding outbursts 

We have developed models to simulate ejected nebulae numeri¬ 
cally and calculate spectra from them at different slit positions. 
The velocity information of the part of the nebula in the slit 
is then converted to a position-velocity plot (PV-plot), which 
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Fig. 3. Ne gative image of Ml-67 made by the Hubble Space 
Telescope (jrosdidier et al. 1998). 



Fig. 4. Explanation of the coordinates and angles used in this 
section. The x and y axes represent the celestial coordinates 
right ascension and declination, the z axis is the radial coordi¬ 
nate, with negative values toward the observer. The thick arrow 
represents the symmetry axis of the structure. Left panel: the 
inclination is changed by rotation around the x axis about the 
angle i. Right panel: the inclined structure is rotated around the 
z axis to change the position angle (f. 


can be compared to observational data. The purpose of these 
models is to give insight in the relation between features in PV- 
plots and the geometry of the nebulae that cause them. Hence 
we have modelled the most common and most likely structures 
for ejected nebulae; spheres, ellipsoids and cones. 

3.1. The geometry of the nebulae 

We start with creating a geometry in the computer. Our struc¬ 
tures are hollow and made up of single dots that represent the 
surface of the nebula. Solid or partially filled structures can be 
made by nesting different shells into each other. For each dot, 
that represents a small volume, the three rectangular coordi¬ 
nates with respect to the star are calculated. 

For spheres we use; 

X - Rneh COS l/r siu 0 1 

y = /;„eb sin sin 0 i (1) 

Z = /;„ebCOS0 j 

where 7?neb is the radius of the nebula, f ranges from 0 to 2n and 
0 from to -i-|, in such way that the dots are evenly spread 
over the spherical surface. 

Ellipsoids are made in the same way. We only consider ax¬ 
ially symmetric ellipsoids. Their geometry is then defined by 
the oblateness m = /?poi/f?eq, with R^\ and /?eq the polar and 
equatorial radius respectively. The x and y coordinates are then 
similar to that of the sphere in Eq. and the z coordinate is 
given by; 

Z = Rneh cos 0 • m (2) 

Bipolar cones are characterised by a semi-opening angle 
(SOA). The nebulae with this geometry can be created using 
the equations; 

.r = z cos i/r tan SOA 1 
y - z sin i/f tan SOA J 


where z varies between -R^eh and H-/?neb and f between 0 and 
In. 

The ellipsoids and cones can be made either elongated to 
simulate jets, or flat to create disks. All structures have axial 
symmetry and are created with the line of sight as their sym¬ 
metry axis (z axis). They are then rotated around a line perpen¬ 
dicular to the line of sight {x axis) to change the inclination 
i and subsequently rotated around the line of sight (which is 
then no longer the symmetry axis) to change the position angle 
ip. This is illustrated in Fig. Q 

3.2. The velocity models 

After choosing the orientation of a structure, we overlay slits of 
finite width and transform the radial coordinates of the underly¬ 
ing points to radial velocities. We assume the velocity-distance 
relation for the nebula to be 



Here, R^eh is the radius of the nebula, vq is the expansion veloc¬ 
ity, r - (x^ +y^ + z^) is the distance to the star and z is the 
radial distance to the star, measured along the line of sight. For 
a we can choose a suitable radius-velocity dependence; a < 0 
for braking after the outburst, a - 0 for a continuous outflow 
with constant velocity, a - 1 for a single, short burst with no 
further interaction and a > 1 for acceleration of the material af¬ 
ter the outburst. For our models of freely expanding outbursts, 
a thin nebular shell with a = 1 is used. Since all our structures 
are made up of single dots, the transformations can be carried 
out point by point. 

For each slit, our models produce one PV-plot, in which the 
radial velocity information under the slit is plotted against the 
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position along the slit. We took all slits parallel to the x axis, so 
that they have a constant y coordinate. The PV-plots thus show 
the radial velocity plotted against the position x. All plots 
are scaled to arbitrary units in position and velocity, because 
we are only interested in the shape of a structure. In actual ob¬ 
servations, the structures may therefore appear stretched out in 
either direction (position or velocity). The result of this whole 
exercise is a ’reference guide’ that we used to determine the 
origin of certain structures in observed long slit spectra]^ 


panding nebula does not hold for Ml-67, is that there is indeed 
a strong interaction between the stellar wind and the ISM and 
that this interaction produces a bow shock around the stellar 
surface. As we will point out in Sect. ^ the bow shock model 
is also qualitatively able to explain the chaotic structure of the 
nebula, which is much more difficult to do with a freely ex¬ 
panding nebula. 

4. The bow shock models 


3.3. Freely expanding outbursts for WR 124 

We tried to fit the different geometries described above and 
combinations of them to the long slit spectra data of Ml-67, 
shown in Fig. The clumpy and chaotic structure of the neb¬ 
ula made it hard to fit. We then also tried to fit parts of the 
nebula. However, we never found a convincing solution. We 
therefore questioned our assumption that the outburst around 
WR 124 is expanding freely and concluded that this may not 
be the case. The Fabry-Perot data of Fig. which we obtained 
afterwards, supported this conclusion. The following features 
in the velocity data of Ml-67 particularly convinced us: 


There is hardly any nebular emission seen at radial veloc¬ 
ities higher than -1-208 km s ', though the star has a ra¬ 
dial velocity of about -1-199 km s *. In case of a continu¬ 
ous outflow by a wind or discrete LBV-like outbursts, one 
would expect to find nebular emission with radial velocities 
of tens, or even hundreds of kilometres per second higher 
than that of the star. 

Emission is seen at radial velocities from about -1-78 
to 4-208 kms^', of which the average is -i-143 kms“*. 
However, the bulk of this emission is found at velocities 
higher than about -1-140 km s“' (See Fig. |). This means 
that the high-velocity part of the nebula shows much more 
emission than the low-velocity pa rt. This was already men¬ 
tioned by Solf & Carsenty ( 1982 ). 

There is an asymmetry in the emission distribution. When 
looking at the different panels in Fig. ^ starting at the high¬ 
est velocity and skipping down, it is clearly seen that the 
emission appears close to the star and becomes increasingly 
broader. When skipping from the lowest velocity panel up¬ 
ward, it is seen that (except for some small patches at lower 
velocities) the emission ’starts’ at about 120-130 km s“* as 
a very broad pattern. The far hemisphere is thus narrower 
than the near hemisphere and the bright emission region 
between -1-140 and -1-208 km s“* looks like a triangle that is 
pointing away from us. 


The explanation that can be found in literature is that the 
nebula as a whole is braked by the ISM, so that the star is dis¬ 
placed from the centre of expansion toward the leading edge 
( |Chu & Treffers 1981; Solf & Carsenty 1982). This means that 
the far and near hemisphere are braked equally, which is in 
contradiction with the asymmetry in the emission distribution. 
Instead, we think that the reason that the model of a freely ex- 


' The simulated long-slit spectra for various geometries and veloc¬ 
ities are available from the first author via e-mail. 


The failure of fitting a freely expanding shell to the velocity 
data of Ml-67, the high radial velocity of WR 124 and the 
global picture from the more detailed Fabry-Perot data show 
that our assumption of the freely expanding shell is incorrect. 
Instead, the high velocity of the star is responsible for the for¬ 
mation of a paraboloid-shaped bow sh ock. The possibi li ty of a 
bow shock was already mentioned by Crosdidier et al. ( 1999 ). 
The bow shock model allows the star to be much closer to the 
ISM, which makes it easier for outbursts to collide with it. We 
will discuss the models for such a bow shock in this section. We 
start with a two-dimensional model and convert it to a rotation- 
ally symmetric three-dimensional model, which can be tilted 
to any wanted orientation, before obtaining the radial velocity 
information. 


4.1. 2D models for continuous mass loss 


The models we use here are published by Canto et al. ( 1996 ). 
They derive analytical expressions to describe wind-wind in¬ 
teractions in general. They also discuss the special case where 
one wind is plane-parallel and the other spherical, as is the 
case when a star with spherical wind is ploughing through a 
homogeneous ISM. In our models we assume the interstellar 
medium to be homogeneous. We also assume that the stellar 
wind is constant over a long period (~ 10^ yr) compared to the 
timescale of LBV-outbursts (usually on the order of 10^ yr), 
except during these outbursts. The wind-ISM interaction is re¬ 
sponsible for the geometry and dynamics of the bow shock, 
while all emission comes from the material that is ejected dur¬ 
ing the outbursts and its collision into the bow shock surface. 
This implies the assumption that the momentum that is carried 
along with the outbursts is lower than the momentum in the 
wind-ISM interaction. 

The distance between the star and the front of the bow 
shock is called the stagnation point distance and can simply 
be derived by momentum equilibrium: 




Mvc 


47TpismVl 


1/2 


( 6 ) 


where Vism is the velocity of the ISM relative to the star (we will 
observe the situation from the rest frame of the star for simplic¬ 
ity). The geometry of the bow shock is expressed in terms of r, 
the distance between the star and a point on the bow shock, as 
a function of the angle 6 that is defined in Fig. ||(a): 

r = ro CSC e^]2>(\-0co^e) 


(7) 
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- 3 - 2-10 1 2 

z coordinate (pc) 


Fig-5- Graphical representation of the bow shock model by 


Canto et al.| ( |1996| ), with the physical parameters for WR 124. 
Upper panel (a): geometry (z versus x coordinate), lower panel 
(b): velocity along the shock surface (v). In both panels the 
horizontal axis represents the symmetry (z) axis of the shock in 
parsecs. 


( [Canto et~ar 1996). The velocity of the material along the sur¬ 
face of the shock is split in a component parallel to the symme¬ 
try axis {z axis, Vz) and a component perpendicular to it (vxy). 
They are given by: 


Vz = 


Mvac sin^ 6 - sin^ 0 


Vxv = 


2M(1 - cos 0) + 47rpisniVismr^ sin^ 0 
Mvoo (0 - sin 0 cos 0) 


2M(1 - cos 0) + 47rpisniVismr^ sin^ 0 

- V" 


/Vz + V^y 


( 8 ) 


(9) 


( 10 ) 


( jCanto et"^ 1996 ), where v is the total velocity along the bow 
shock surface. In Fig. |] the geometry and the total velocity 
along the surface are displayed. 


4.2. 3D models for continuous mass loss 


The model of Canto et al. (1996) is two-dimensional and can 


( 11 ) 


easily be converted into a 3D-model by 

X - rsin0cosi/f) 
y - rsinffsini/^ > 
z = rcosfl j 

where 0 is defined in Fig. ||(a) and ij/ takes evenly spread val¬ 
ues between 0 and 2ji around the symmetry axis. Thus, the 3D 
surface is obtained by rotating Fig. ^(a) around the horizontal 
axis. The velocity components are given by 


= Vx 


►•xy COS \f/ 

- Vxy sin i/f 


( 12 ) 


and by Eq. (||). 

The three-dimensional bow shock model is now complete 
and oriented in such a way that the observer is looking into the 
hollow shock from the rear. We can now rotate the shock to 
any orientation we like, about the angles i and (f as indicated 
in Fig. The velocity components are rotated the same way, 
and we can use the celestial coordinates x and y and the radial 
velocity Vz to compare the predicted radial velocity maps to the 
observations. 


Fig. I displays the Right Ascension against the radial ve¬ 
locity Vz of a number of bow shocks with different orienta¬ 
tions. From these plots, two effects of the orientation of the 
bow shock can easily be seen. Firstly, when looking at the dif¬ 
ferent columns from the left to the right, the inclination of the 
bow shock increases. The most important result of this is that 
the maximum radial velocity observed shifts from 0 km s * 
(which is the radial velocity of the star) for i = 0° to more than 
30 km s ' for i = 50°. This means that the velocity difference 
between the maximum radial velocity where nebular emission 
is seen and the radial velocity of the star puts an upper limit 
to the inclination. As we have seen in Sect. 3.3, in the case of 
Ml-67 this difference is approximately 10 km s^*, which gives 
a strong indication for a small inclination of the main axis of 
the bow shock. Secondly, when one looks through the different 
rows from the top to the bottom, the position angle changes. 
The result is that the axial symmetry around RA = 0 disap¬ 
pears. We will later see that for the case of Ml-67, only about 
the upper 100 km s * in radial velocity are observed, so that 
this effect manifests itself as a steep drop in radial velocity at 
one side of the maximum and a slower drop at the other side. 
The side of the steep drop indicates the direction where the top 
of the bow shock points at. For example, at the lower right plot 
of Fig. ^ the steep side is at the right, which means that the 
top of the bow shock is pointing to this side, so that the star 
is moving, with respect to its local ISM, toward higher Right 
Ascension. 


4.3. Orientation of the bow shock 

It is possible to roughly say something about the orientation of 
a bow shock of a given star with respect to the line of sight 
(LOS). This way, we can limit the total number of possible 
orientations of the bow shock in Ml-67 appreciably. In order 
to do so, we need to know the velocities of both the star and 
the ISM surrounding the star. The orientation of the main axis 
of the bow shock is simply the orientation of the vector that 
is the difference of the spatial velocity of the star and that of 
the local ISM. The spatial velocity of the star can be derived 
from the measured radial velocity and proper motion. For the 
ISM we will assume that it follows the laws of galactic rotation 
and derive its velocity relative to the sun. After subtracting the 
ISM velocity from that of the star, we can calculate the total 
velocity difference and the inclination and position angle of the 
main axis of the bow shock. We will derive these quantities 
here for the general case of a star with galactic longitude and 
latitude I and b respectively and distance d from the sun, and 
then apply this to Ml-67 in Sect. 
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Fig. 6. Model output for the 3D bow shock models. For all models Vism = 200 km s“* and the star is placed at a Right Ascension 
of 0" and a radial velocity of 0 km s“*. The different rows have position angles ((p) of 0,45 and 90°, the columns have inclinations 
(/) of 0, 10, 20, 30, 40, and 50°. For each panel, the horizontal axis shows the Right Ascension (RA) in arcseconds, the vertical 
axis displays the radial velocity in km s '. See the text in Sect. 4.2 for more explanation and interpretation. 



Fig. 7. The Milky Way as seen from the NGP, with the defini¬ 
tions of the quantities used in this section. 


In the observations of the radial velocity and proper motion 
of the star the motion of the sun is already taken into account, 
so that we indeed need to correct for the motion of the sun when 
calculating the velocity of the ISM. First, we will correct for the 
motion of the local standard of rest (LSR), that moves with a 
velocity of vq - 220 km s^* around the galactic centre (GC) 
and subsequently the peculiar motion of the sun with respect 
to the LSR will be taken into account. This peculiar velocity is 
calculated from Hipparcos ob servations for d ifferent classes of 
stars in the galactic plane by |Vlignard ( 200C ). We will use the 
averages of the values given in that article: 


uq - 10.0+1.3 km s * I 

v© = 13.9+ 3.7 km s^' i (13) 

W 0 = 7.4 + 2.6 km s“* j 


for the direction toward the GC, the direction of galactic rota¬ 
tion and toward the north galactic pole (NGP) respectively. 

When looking along the line of sight (LOS) defined by the 
galactic coordinates / and b, we can calculate for a point at any 
distance from the sun d, the distance from that point to the GC, 
which we will call R (see Fig. 

R^[d^sir?l + {RQ-dcoslf)^^^ (14) 
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where Rq is the standard value for the distance from the sun 
to the GC; Rq = 8.5 kpc. We can then calculate the angle y, 
defined in Fig. by 


y = arctan 


Rq - d cos I 
t/sin/ 


(15) 


The circular velocity for a distance R with 3 kpc < R < Rq is 
given by |Burton| (|1988|): 


VQ 


L0074|A)' 


0.0382 


0.00698 


(16) 


The components of the ISM velocity in the radial direction and 
the direction of galactic longitude and latitude, are then given 
by: 

Vism,r = Vc COS(y - 1) COS b \ 

Vism,i = Vc sin(y -/) COS Z? I (17) 

t'ism.b “ t^c sin b j 


In these expressions, we have neglected the velocity component 
of the ISM perpendicular to the plane of the galaxy. 

The proper motion of the star can be expressed in terms of 
galactic coordinates (jU/ and fih). The proper motion gives rise to 
spatial velocity components that are dependent of the distance. 
We thus have: 


V»,r 

V.,1 =;«;■£/ i (18) 

V.,b = Rb-dj 

We can now calculate the velocity difference between the 
star and the local ISM in the three components, where we also 
correct for the motion of the sun: 

Vr = v»j - (vism.r “ VQ sin I - uq COS I - VQ sin 1) 1 

Vl - V,,l - (Vism.l - Vo cos 1 + Mq sin 1 - Vq cos /) > (19) 

Vb = V,,b - (Vism,b + Wq) j 


The total velocity of the star with respect to the local ISM 
is what we called Vism in Sect. 4.1: 


= ( 


V, + V, + v. 


( 20 ) 


Furthermore, we can also find an expression for the inclination 
i and position angle ip: 


arccos 





- arctan 



+ e 


( 21 ) 

( 22 ) 


where e denotes the angle between the line of constant declina¬ 
tion and the line of constant galactic latitude at the position of 
the star. 

The uncertainties in the distance and in the proper motion 
of the star introduce a range of possible values for vism, i and ip, 
which is a constraint that we can use to reduce the total num¬ 
ber of possible orientations drastically. This makes it easier to 
compare our models to the observations, as we will do in the 
next section. 


E 



-60 -40 -20 0 20 40 60 80 

Right Ascension (arcsec) 


G 


-60 -40 -20 0 20 40 60 

Declination (arcsec) 



Fig. 8. Negative images of the emission distribution as inferred 
from the Fabry-Perot images in Fig. Upper panel (a): radial 
velocity against Right Ascension, Lower Panel (b): radial ve¬ 
locity against Declination. The process of creating this image 
is described in the text. The disturbances in the edges of the im¬ 
age are artifacts caused by manipulating the image. The sharp 
cutoff in the lower image at abou t -20 "is caused by the seeing 
during the observations (see Sect. 2.2). 


5. Comparison of bow shock models for 
continuous mass loss with observations 


In this section we will use the bow shock model discussed in 
Sect. ^ to explain the distribution of emission as observed in 
the Fabry-Perot images. In order to do this properly, we trans¬ 
formed the observational data in Fig. to an image that dis¬ 
plays the emission distribution in the Right Ascension - radial 
velocity plane. This was done by stacking the 30 different FP 
images on top of each other to create a 3D body, with two spa¬ 
tial dimensions and one radial velocity dimension. Then we 
summed all layers with constant declination to get the image 
that is shown in Fig. ||(a). Fig. ^b) shows the same image, but 
here the Right Ascension is replaced by the declination. These 
are the images we will compare our bow shock models with 
different parameters to. 


We apply the properties of WR 124 to the results that we 
derived in Sect. 4.3. The star has galactic coordinates / = 50.2°, 
b = 4-3.31°. At this position, the angle e has the value of 62.7°. 
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Fig. 9. Theoretical properties of the orientation of the bow 
shock around WR 124 as a function of the distance from the 
sun: Upper panel: Total velocity of WR 124 relative to the ISM. 
Middle panel: Inclination (/) of the main axis of the bow shock. 
Lower panel: Position Angle {(p) of the main axis of the bow 
shock. The solid line gives the most probable value, the dashed 
lines display the uncertainties that are a result of the uncertainty 
in the proper motion of the star. 


The resemblance between the models and the data is re¬ 
markable. It is clear that the best model can be found some¬ 
where in the hrst row and the third or fourth column of Fig. 
From that, it can be inferred that the best values for the param¬ 
eters are Vism 180 kms^*, i 20° and (f -185°, so that 
the star is moving south, under a very small inclination with 
respect to the line of sight. The resemblance between the data 
and the model indicates that the values for these three param¬ 
eters ht well, not that this distance is better than other values. 
The distance was only used to obtain likely values for these pa¬ 
rameters. As can be seen from Fig. there are more distances 
that can explain this combination of the three parameters. We 
conclude that the bow shock model with the parameters given 
above is successful in explaining the overall features of the two 
velocity versus position (a or 6) plots of Fig. || and Fig. |^. 

There are of course also differences between the model out¬ 
put and the observations. Part of these differences can be at¬ 
tributed to the fact that we assume a continuous mass loss from 
the star, which is in reality unlikely to be the case. In fact, in 
Sect. 1^ we will argue that the mass loss occurred in outbursts. 
Also, we assume the ISM to be homogeneous. Inhomogeneities 
in the ISM will cause a less smooth bow shock surface, as in¬ 
deed seen in the observations. The largest disagreement be¬ 
tween model and observations is of course seen in the lower 
panel of Fig. |^. This is due to the lack of observational data 
for declinations less than about -20"(see Sect. 2.2 and Fig. ||). 


The proper motion of WR 124 was measured by Hipparcos and 
can be converted to galactic coordinates: 


Pi = -6.1 + 2.0 mas yr ’ 
Pb - -3.0 + 2.0 mas yr“' 


(23) 


We can then use Eqs. (|2^, ( ^1| ) and to calculate the 
velocity difference between the star and the ISM (vism), as well 
as the inclination (/) and position angle {ip) of the main axis 
of the bow shock (as dehned in Fig. as a function of the 
distance {d) of WR 124 from the sun. The results are plotted 
in Fig. ^ where the most probable value, the upper and lower 
limit of the proper motion of Ml-67 are depicted in the solid 
and dashed lines respectively. Because of the great distance to 
the star, the uncertainties in the proper motion are quite large, 
so that the resulting uncertainties in the velocity components 
perpendicular to the LOS are also large and increase with the 
distance. 


Next, we apply the results from Sect. 4.1 For d - 6.5 kpc, 
Eq. gives Vism ~ 174 km s * and we hnd from Eq. (H) that 
ro = 1.3 pc, where we used M - 2.45 x 10“^ Mq yr“', Voo = 
710 km s^* ( Nugis & Earner^ 200C) andpism = 1 niyi cm 

We can now build the 3D models as described in Sect. |4.2| 
and rotate them according to the angles found in Eig. y. 
Because of the large uncertainty in the proper motion, we cal¬ 
culate three different models for each distance d, with the three 
different values found at each line in Eig. ^ for Vism, i and tp 
and plot a contour of the model output over Eig. |8[ The result 
is displayed in Eig. for distances of 6.0, 6.5, 7.0, 7.5, 8.0, 
and 8.5 kpc. 


6. The effect of outbursts on the bow shock 

So far, we have looked at a bow shock that is formed due to a 
continuous wind. However, as we have already noted in Sect. [| 
and 2.3, the distinct arcs that can be seen in the HST image are 
likely to originate from discrete outbursts. Since these kinds 
of outbursts only take place after the star has moved from the 
main sequence and since the O-star that preceded WR 124 must 
already have had a strong wind, these outbursts must have oc¬ 
curred inside the bow shock. In this section we will describe 
what the effects of a short outburst on the bow shock is. 

In the hrst phase the outburst will expand freely and both 
the geometrical and velocity structure are undisturbed, apart 
from the extra radial velocity due to the velocity of the star. 
Typical timescales for this phase are derived from division of 
the distance ro of a few parsec by a typical LBV outburst ve¬ 
locity of 100 kms“*, which gives a few times 10"^ yr. After 
this time, the outburst will impact onto the bow shock surface, 
which will affect the observed nebula in different ways. Eirstly, 
the outburst will be braked or possibly even be halted by the 
surface of the bow shock. This means that the radial velocity 
will change drastically and that the nebula will brighten. Even 
in the case of a spherical outburst, the different parts of the out¬ 
burst will reach the bow shock at different times, so that these 
effects of braking and brightening will propagate through the 
nebula, starting at its top. Secondly, the gas of the outburst will 
be dragged along the surface of the bow shock and eventually 
adopt its velocity. This means again a difference in the veloc¬ 
ity pattern and a distortion in the (thus far possibly symmetric) 
geometry of the outburst. 
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Fig. 10. Model output for the 3D bow shock models for WR 124. Upper panel (a): Right Ascension (") versus radial velocity 
(km s“'), lower panel (b): Declination (") versus radial velocity (km s“'). The grey-scaled images in the background are smoothed 
versions of Fig. The thin lines are contours of the data. The model output is plotted as the thick curve, that encloses the actual 
bow shock in this projection and that would look like the output in Fig. ^ if plotted completely. The different columns have values 
for the total velocity, inclination and the position angle for the distances of 6.0, 6.5, 7.0, 7.5, 8.0 and 8.5 kpc. The rows are the 
three different values for these quantities for the lower limit, the best value and the upper limit for the proper motion respectively, 
that each row displays one line of each graph in Fig. The dashed lines in each plot are at 0" and 199 km s ', the position 


so 


and radial velocity of WR 124. Note the lack of observational data at declinations lower than -20"(See Sect. 2.2). 


Because the magnitude of the effect described above will 
depend strongly on the relative momenta of the matter in both 
the bow shock and the outburst, we divide the outbursts in three 
simple categories: 

1. Spherical outbursts where the momentum of the outburst is 
smaller than that of the wind-ISM shock 


2. Spherical outbursts where the momentum of the outburst is 
larger than that of the wind-ISM shock 

3. Non-spherical outbursts that may be a combination of the 
two cases above: in some directions ^outburst < Fbowshock. in 

others p outburst ^ Pbowshock- 
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In the first case, the geometry and dynamics of the bow 
shock will not be changed too much by the impact of the out¬ 
burst. Those of the outburst, however, will change apprecia¬ 
bly. The freely expanding gas will impact onto the bow shock, 
brighten and be dragged along the surface layer of the bow 
shock. The outburst will adopt the geometry and dynamics of 
the bow shock quickly, so that it is eventually only recognis¬ 
able as a bright, distorted ellipse in the diffuse, weak, if at all 
visible, background of the gas that resides in the bow shock. 

In the second case, the steady bow shock is distorted by the 
impact. The outburst will be brightened and slowed down at 
the moment of impact, but will continue its motion for a while. 
Eventually, the outburst will be halted by the ongoing ram- 
pressure of the ISM, but this will happen at a larger distance 
from the star than where it hrst hit the bow shock. This situa¬ 
tion can be described in the same way as the situation before 
the outburst, only with a stronger wind (with lower velocity but 
a much higher mass loss). This will result in a larger stagnation 
point distance ro, as given by Eq. (||). The outburst will change 
the shape of the bow shock temporarily, because the different 
parts of the outburst will impact at different times. The distor¬ 
tion will thus propagate through the bow shock, starting at its 
top. Eventually, this results in a bow shock with the same ge¬ 
ometry, but a different size and different dynamics. The size 
simply scales with ro, but the velocities in Eq. (Q) and Eq. 
depend on the mass loss and terminal wind velocity in a more 
complex way. After the outburst, the momentum of the quiet 
wind will drop back to its old value, because the outburst itself 
lasts only very short compared to the evolutionary timescale of 
the star in this phase. Because of this, the old value of ro may 
be reinstalled after some time. 

In the third case a non-spherical outburst occurs. Many 
LBV outbursts are known to be bipolar, the best known exam¬ 
ple being rj Carinae. Bipolarity of an LBV outburst can result 
from binarity, or from fast rotation. Since an interacting spheri¬ 
cal outburst as described above can at best result in an elliptical 
ring, the structures in this case are already more complex. This 
complexity will increase when the impact occurs. In the case 
of a non-spherical outburst, it is likely that the outburst will not 
hit the bow shock at its top first, but somewhere on the side, 
where the surface velocity of the shock is much higher. Thus 
the resulting drag will be stronger and will decrease symmetry 
even further. 

In all cases, part of the outburst will move toward the rear 
end of the bow shock. Because this end is open, it will never 
interact with the bow shock surface and keep expanding freely. 
This part of the outburst can therefore be used to derive a dy¬ 
namical timescale directly. In the velocity data of Ml-67, we 
hnd a freely expanding structure that is labelled as A in Fig. |^. 
From comparison with Fig. it is clear that this structure is 
not located on the bow shock surface. We fitted the structure A 
to an ellipse, assuming that WR 124 is at the centre of expan¬ 
sion. From the ht, we found Vexp ~ 150+15 km s“* and r =» 40" 
(oc 1.3 pc at £/ = 6.5 kpc), which gives a dynamical timescale 
of 1250 + 125 yr kpc“* or about 8.2 ± 0.8 x 10^ yr for for the 
estimated distance of 6.5 kpc for WR 124. This timescale is 
typical for LBV outbursts and therefore a strong hint that the 
nebula Ml-67 may indeed be the result of LBV-like eruptions. 



60 40 20 0 20 40 60 

Riyhl Ayceiisioii (arc«ec} 


Fig. 11. The two data sets used for this research. Upper left 
panel (a): HST image, upper right panel (b): Fabry-Perot data 
projected in the declination-radial velocity plane, lower left 
panel (c): Fabry-Perot data projected in the right ascension- 
radial velocity plane. In the images, parts of the nebula that are 
discussed in the text are encircled and labelled, for easier refer¬ 
ence. The reader can also easily compare the positions of each 
structure in the different images. Note that (roughly) the part of 
the upper right image where no data is available is left blank. 


It is very unlikely that the structure A is the only result of the 
outburst that caused it, which we will refer to as outburst A. 
The structure is clearly expelled toward us, and one would ex¬ 
pect to find a counterpart of this structure that was caused by 
the same outburst, but directed away from us. This counterpart 
should be clearly visible, as long as it hasn’t reached the bow 
shock surface yet, because it should have a much higher radial 
velocity than the star (up to +345 km s“*). Since we do not 
see this counterpart, it must have collided with the bow shock 
surface already. Given the expansion velocity derived above 
(150 km s“') and the stagnation point distance of ro ~ 1.3 pc 
(Sect. it would take the part of outburst A that is directed 
away from us about 8.5 x 10^ yr to reach the bow shock sur¬ 
face. Within the accuracy of the ht, this time agrees with the 
htted dynamical timescale of 8.2 + 0.8 x 10^ yr. This means 
that the counterpart of the freely expanding structure A might 
have collided with the bow shock just a short time ago. It is 
therefore very likely that the structure A originates from the 
most recent outburst that occurred on WR 124. 

The interactions between parts of an outburst and the steady 
bow shock can cause a very chaotically looking nebula, just as 
is the case for Ml-67. In particular, it will generate elliptical 
structures of which many can be seen in Fig. ^ Another result 
of the impact of an outburst onto the bow shock, is that it is no 
longer possible to simply fit an undisturbed, freely expanding 
nebula to this part and derive its dynamical timescale. However, 
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we may still be able to estimate the dynamical timescales of 
parts of the outburst that have been interacting. 

The majority of the emission of Ml-67 seems to originate 
from the surface of the bow shock in the position-velocity di¬ 
agram of Fig. 1^. Consider the structure that we labelled B in 
Fig. 0. From the Fabry-Perot images, we derive its velocity, 
as shown in Fig. |^c). This figure, compared to the hrst row, 
fourth column image of Fig. |^a), shows that the structure lies 
on the bow shock surface. We do not know the exact path that 
it travelled to get to its current position, but we do know that 
this path must be in the plane that is defined by the main axis 
of the bow shock and the current position of this structure B. 
We can therefore characterise the path by a single angle 9, that 
is the angle between the line that joins the front of the bow 
shock and the star and the direction in which the structure was 
expelled, as depicted in Fig. |^a). After the outburst B, that 
caused the structure B, occurred, the ejecta hrst crossed the dis¬ 
tance between the star and the bow shock surface in a straight 
line and at a constant velocity. After having reached the bow 
shock surface, the ejecta collided with it and were eventually 
dragged away along the bow shock surface toward the rear, un¬ 
til they reached their current position. The angle 6 can have 
values ranging from 0° (expelled exactly toward the front of 
the bow shock) to about 70°, the angle at which the structure is 
observed now, which would mean that it moved from the star 
straight to its current position. If we assume that the structure B 
was expelled at the same outburst as the freely expanding struc¬ 
ture A discussed above, it should have crossed the distance to 
the bow shock surface with a velocity of Vexp = 150 km s”'. 
The velocity on the bow shock, along its surface, is given by 
Eq. (p^. If we then integrate over the path, we hnd a dynam¬ 
ical timescale as a function of 9. The results are displayed in 
Fig.[T^b). From this, we find that the dynamical timescale for 
structure B must be greater than 1.2 x 10"^ yr, for the case where 
this structure has moved directly from the star to its current po¬ 
sition, and that it is very likely that the dynamical timescale is 
in the order of 2 x 10“* yr (for 9 45°). For this 9, the timescale 

would even increase to about 4x10"^ yr, if the expansion ve¬ 
locity would be 50 kms^'. For very small 9, the dynamical 
timescale goes to inhnity, because the velocity alo^ the bow 
shock surface drops to zero near the front (See Fig. ^(b)). 

One could debate whether the ejecta can survive the harsh 
stellar environment for a few times 10^ yr. Indeed, we see that 
blobs of gas around young, hot stars are photoionised by the 
strong radiation held of the star, after which the ionised gas 
is blown away by the stellar wind. Being a Wolf-Rayet star, 
WR 124 clearly has such a strong radiation held. However, 
the gas cannot been blown away by the wind, since there is 
an equilibrium on the surface of the bow shock between the 
stellar wind and the ram pressure of the ISM. The gas therefore 
cannot escape, other than in the direction toward the back of 
the bow shock, as we described. In fact, the emission from Ml- 
67 that we see in Fig. ^ is in Ha, which indicates that the gas 
is already ionised. Since we only consider the dynamics of the 
gas here, it is not important what exactly happens to the gas, as 
long as it does not affect its global motion. 

For structure C in Fig. |^, roughly the same scenario holds 
as for structure B. It is very likely to be located on the bow 


L5 
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Fig. 12. Upper panel (a): Geometry plot, showing the possible 
paths from the star (asterisk at the origin) to an observed struc¬ 
ture of gas on the bow shock surface (asterisk at left), char¬ 
acterised by the angle 9. Lower panel (b): Travel time for an 
ejected structure to reach the observed position as a function of 
9. The dashed line is the time needed to cross the distance from 
the star to the bow shock. We assumed an expansion velocity of 
v'exp = 150 km s *. The dash-dotted line is the time needed to 
travel the second part of the path, along the bow shock, where 
the velocity is given by Eq. (|l0|). The solid line is the total travel 
time. 


shock surface. However, this structure it spread out much more 
and therefore it is hard to define its exact location. Eurthermore, 
it is positioned at the southern cutoff limit of our velocity data. 
Eor these reasons, we do not discuss its dynamical properties 
here. We labelled it in Eig. [l^ for completeness, because it is 
another feature that is clearly visible in the different datasets. 

Although there are still large uncertainties in the rough 
calculations presented here, it seems very likely that the dy¬ 
namical timescales of the freely expanding structure A and the 
shocked structure B do not match. The explanation that both 
structures originate from the same bipolar outburst, with an ex¬ 
pansion velocity that is different in different directions seems 
unlikely, since the directions to these two structures as seen 
from the star do not differ too much. This means that the two 
features are likely to originate from different outburst events. 
This result strongly suggests a multiple-outburst scenario, with 
typical outburst timescales of a few times 10"^ yr. This scenario 
favours an LBV phase rather than an RSG phase that is respon¬ 
sible for the nebula Ml-67. 
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7. Summary and conclusions 

We tried to fit different freely expanding models to the nebula 
Ml-67 and find that Ml-67 is not expanding freely, but that 
the outbursts that formed the nebula have been interacting with 
the ISM. Evidence for this is given by the fact that the nebula 
has a radial velocity that is lower than the radial velocity of 
the star, in particular the fact that there is no nebular emission 
with a significantly higher radial velocity than that of the star. 
The asymmetry in the amount of emission between the high- 
velocity and low-velocity sides of the nebula and the asymme¬ 
try in the velocity distribution of the nebula seem to confirm 
the fact that there are interactions between the outbursts and 
the ISM. 

We find that the central star WR 124 has a velocity of 
about 180 km s * with respect to the surrounding ISM, mov¬ 
ing roughly away from us, which causes a paraboloid-like bow 
shock around the star. Because the star is moving away from 
us, we are looking into the hollow bow shock from the rear. 
Using the simple assumption of momentum conservation, we 
can model the geometry of the bow shock and the velocity 
of the gas that is moving along its surface. This way we can 
roughly derive the orientation of the bow shock, which is de¬ 
termined by the spatial velocities of the star and that of the 
ISM surrounding it. The inclination of the main axis of the bow 
shock with respect to the line of sight is about 20°, where the 
arrow that points from the star to the front of the bow shock 
is mainly pointed away from us and slightly to the south. This 
small inclination is responsible for the fact that we cannot dis¬ 
tinguish the bow shock in the HST image in Fig. ^ and that 
we see the star projected in the centre of the nebula. The star 
sits inside this bow shock, at about 1.3 pc from its front. As 
input for our bow shock models, we used the observed prop¬ 
erties of the stellar wind (M, Voo) and assumed an ISM density 
of 1 mn cm“^. These values obviously have quite some uncer¬ 
tainty in them. We therefore once more explicitly mention the 
rough nature of these calculations. The assumptions of a con¬ 
tinuous wind and a homogeneous ISM for our models are very 
likely to be incorrect. In fact, an ISM with inhomogeneities 
could explain why the observed bow shock does not have a 
smooth surface. 

In general, a run-away star with a wind will form a bow 
shock. The LBV eruptions during the late stages of the evo¬ 
lution of massive stars will partly collide with the bow shock 
surface, brighten up and be dragged away by the ram pressure 
of the ISM, along the surface of the bow shock. In this way, 
irregular, or even chaotic nebulae with arcs and rings can be 
formed by intrinsically nicely behaving, possibly even spher¬ 
ical, discrete outbursts. The interactions with the ISM make 
the derivation of dynamical timescales for the outbursts much 
more difficult than in the case of low-velocity stars, which blow 
voids of tens of parsecs in the interstellar medium during their 
O-star phase, so that outbursts occurring in a later stage will in¬ 
deed expand freely. However, also in the case of a bow shock, a 
part of each outburst will move toward the rear end of this bow 
shock and remain expanding freely forever. This part might still 
be fitted in the straightforward way. 


For Ml-67 we were able to fit the freely expanding part of 
the outburst A in Fig. to r 40" and Vexp ~ 150 km s“', 
which gives the rough dynamical timescale of 8 x 10^ yr, as¬ 
suming a distance of 6.5 kpc. This dynamical timescale is of the 
same order as the time it would take this outburst to cross the 
distance of 1.3 pc from the star to the front of the bow shock, 
which means that the far part of this particular outburst might 
just have collided with the bow shock. In fact, the collision of 
this part of the nebula is required in order to explain the lack 
of nebular emission with significantly higher radial velocities 
than that of the star. 

Part of the nebular emission is found to be on the bow shock 
surface. In order to get there, this matter could have followed 
different paths, because we do not know in which direction the 
material has been ejected from the star. We can calculate the 
amount of time for this matter to have reached its observed 
location, as a function of the path. We did so for the arc that 
is labelled B in Fig. assuming that the expansion velocity 
was the same as for the freely expanding structure A, Vexp = 
150 km s *. We then find that the age of structure B must be 
greater than 1.2 x 10^ yr, probably in the order of 2 x 10^ yr, 
so that this arc most likely originates from a different outburst 
than structure A. 

The multiple outbursts, the expansion velocity of about 
150 km s ' and the dynamical timescales in the order of lO"' yr 
are all typical for LBVs. We therefore conclude that an LBV 
phase preceded the current Wolf-Rayet phase. The outbursts 
that occurred during this LBV phase seem to be the most likely 
explanation for the creation of the Wolf-Rayet ring nebula Ml- 
67 around WR 124. 

Acknowledgements. We thank Thierry le Fevre for his contribution to 
the freely expanding wind models, Antonella Nota for the long-slit 
spectra and Yves Grosdidier for the Fabry-Perot images. 

References 

Bertola, F. 1964, PASP, 76, 241-H 

Burton, W. B. 1988, in Galactic and Extragalactic Radio 
Astronomy, 295-358 

Canto, J., Raga, A. C., & Wilkin, F. P. 1996, ApJ, 469, 729 -h 
C hu, Y.-H. & Treffers, R. R. 1981, ApJ, 249, 586 
Grosdidier, Y, Moffat, A. F. J., Joncas, G., & Acker, A.. 1998, 
ApJ, 506, LI27 

Grosdidier, Y, Moffat, A. F. J., Joncas, G., & Acker, A. 1999, 
in ASP Conf. Ser. 168: New Perspectives on the Interstellar 
Medium, 453-1- 

Lamers, H. J. G. L. M. & Cassinelli, J. P. 1999, Introduction 
to stellar winds (Introduction to stellar winds / Henny 
J.G.L.M. Lamers and Joseph P. Cassinelli. Cambridge ; 
New York : Cambridge University Press, 1999. ISBN 
0521593980) 

Lamers, H. J. G. L. M., Nota, A., Panagia, N., Smith, L. J., & 
Langer, N. 2001, ApJ, 551, 764 
Langer, N., Hamann, W.-R., Lennon, M., et al. 1994, A&A, 
290, 819 

Merrill, P. W. 1938, PASP, 50, 350 h- 
Mignard, F. 2000, A&A, 354, 522 


Van der Sluys &. Lamers: The dynamics of the nebula Ml-67 around the star WR 124 


Minkowski, R. 1946, PASP, 58, 305 + 

Nota, A. & Clampin, M. 1997, in ASP Conf. Ser. 120: 

Luminous Blue Variables: Massive Stars in Transition, 303+ 
Nugis, T. & Lamers, H. J. G. L. M. 2000, A&A, 360, 227 
Pismis, P & Recillas-Cruz, E. 1979, Revista Mexicana de 
Astronomia y Astrofisica, 4, 271 
Sirianni, M., Nota, A., Pasquali, A., & Clampin, M. 1998, 
A&A, 335, 1029 

Solf, J. & Carsenty, U. 1982, A&A, 116, 54 
Stothers, R. B. & Chin, C. 1993, ApJ, 408, L85 
—. 1996, ApJ, 468, 842+ 



